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Abstract 

We calculate the precession of a gyroscope at rest in a Bondi space- 
time. It is shown that, far from the source, the leading term in the 
rate of precession of the gyroscope is simply expressed through the 
news function of the system, and vanishes if and only if there is no 
news. Rough estimates are presented, illustrating the order of mag- 
nitude of the expected effect for different scenarios. It is also shown 

from the next order term ( — ^) that non-radiative (but time depen- 
dent) spacetimes will produce a gyroscope precession of that order, 
providing thereby "observational" evidence for the violation of the 
Huygens's principle. 
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1 Introduction 



The theoretical description and the experimental observation of gravitational 
radiation are among the most relevant challenges confronting general rela- 
tivity. 

A great deal of work has been done so far in order to provide a consistent 
framework for the study of such phenomenon. Also,since Weber's pioneering 
work [|l| important collaboration efforts have been carried on, and are now 
under consideration, to put in evidence gravitational waves (see [0, P[, PJ,[§], 
and references therein). 

It is the purpose of this work to evaluate the influence of gravitational ra- 
diation on a gyroscope. This idea is not new, in fact some years ago, Chaboyer 
and Henriksen || put forward the possibility of detecting gravitational radia- 
tion by means of an orbital laser gyroscope. In such experiment the presence 
of gravitational radiation is brought out by the differential effect that radia- 
tion has on the paths of the photons in the rotating frame of reference. 

In this work we shall calculate the rate of precession of a gyroscope in the 
field of gravitational radiation. To do so we shall use the Bondi's formalism 
[0] which has, among other things, the virtue of providing a clear and precise 
criterion for the existence of gravitational radiation (see also Q). Namely, if 
the news function is zero over a time interval, then there is no radiation during 
that interval. Also, the present approach has the advantage of providing a 
very simple expression linking an "observable" (at least in principle) quantity 
(the rate of precession of a gyroscope) with the emission rate of gravitational 
radiation. 

The formalism has as its main drawback || the fact that it is based on a 
series expansion which could not give closed solutions and which raises unan- 
swered questions about convergence and appropriateness of the expansion. 

However since we shall assume the gyroscope to be very far from the 
source, we shall use in our calculations only the leading terms in the expan- 
sion of metric functions. Furthermore, since the source is assumed to radiate 
during a finite interval, then no problem of convergence appears . 

We shall see that the leading term of the rate of precession of the gyro- 
scope (Q) is expressed through the news function in such a way that it will 
vanish if and only if there is no news (no radiation). 

In the special case of the quadrupole radiation (in the linear approxi- 
mation), the rate of precession may be expressed through the third time 
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derivative of the quadrupole moment, or alternatively, through the rate of 
loss of the mass function. ^ 

Next we present the order (— ) of Q. As we shall see, it contains terms 

with news, together with a time dependent term not involving news. This 
last term represents the class of non-radiative motions discussed by Bondi J7| 
and may be thought to correspond to the tail of the wave, appearing after 
the radiation process 0. The obtained expression allows for "measuring" 
(in a gedanken experiment, at least) the wave-tail field. This in turn implies 
that observing the gyroscope, for a period of time from an initial static situ- 
ation until after the vanishing of the news, should allow for an unambiguous 
identification of a gravitational radiation process. 

In the next section we briefly present the Bondi's formalism. The ex- 
pression for the rate of precession of the gyroscope in the Bondi metric is 
calculated in sections 3 and 4, and estimates for the rate of precession in dif- 
ferent scenarios are presented in section 5. Finally the results are discussed 
in the last section. 

2 The Bondi's Formalism 

The general form of an axially symmetric asymptotically flat metric given by 
Bondi is 

ds 2 = (T e 2l3 -U 2 r 2 e 2 ^du 2 + 2e 2(3 dudr 

+ 2Ur 2 e^dud9 - r 2 (e 2l d9 2 + e" 27 sin 2 9d<j) 2 ) (1) 

where V, /3, U and 7 are functions of u, r and 9. 

We number the coordinates x ' 1,2 ' 3 = u, r, 8, respectively, u is a timelike 
coordinate such that u = constant defines a null surface. In flat spacetime 
this surface coincides with the null light cone open to the future, r is a 
null coordinate (g rr = 0) and 9 and 4> are two angle coordinates (see [0] for 
details). 

Regularity conditions in the neighborhood of the polar axis (sin# = 0), 
implies that as sin 9— > 

V, p,U/ sin #,7/ sin 2 9 (2) 
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each equals a function of cos 9 regular on the polar axis. 

The four metric functions are assumed to be expanded in series of 1/r, 
then using field equations Bondi gets 

7 = c(u, 9)r~ l + (C(u, 9) - ^c 3 ) r' 3 + ... (3) 
U = - (c e + 2c cot 9) r~ 2 + ... (4) 

V = r-2M{u,9) 

- (^N e + N cot 9 -c 2 9 -4cc e cot 9 -^c 2 (1 + 8 cot 2 9)J r' 1 + ... (5) 

(3 = -\c 2 r- 2 + ... (6) 
where letters as subscripts denote derivatives, and 

AC U = 2c 2 c u + 2cM + N cot 9 - N e (7) 

The three functions c, M and N are further related by the supplementary 
conditions 

M u = -c 2 u + - (c ee + 3c e cot 9 - 2c) u (8) 

- 3N U = M e + 3cc u9 + Acc u cot 9 + c u c e (9) 

In the static case M equals the mass of the system whereas iV and C are 
closely related to the dipole and quadrupole moment respectively. 
Next, Bondi defines the mass m(u) of the system as 

m( u ) = r M{u, 9) sin BdO (10) 
2 jo 

which by virtue of (H) and @) yields 

m u = -- r clsm9d9 (11) 



2 Jo 

Let us now recall the main conclusions emerging from the Bondi's ap- 
proach. 
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1. If 7, M and N are known for some u = a(constant) and c u (the news 
function) is known for all u in the interval a < u < b, then the system is 
fully determined in that interval. In other words, whatever happens at 
the source, leading to changes in the field, it can only do so by affecting 
c u and viceversa. At the light of this comment the relationship between 
news function and the occurrence of radiation becomes clear. 

2. As it follows from (|TT|) , the mass of a system is constant if and only if 
there are no news. 

In the next section we calculate the rate of precession of a gyroscope at 
rest in the frame of (|1]) 



3 The gyroscopic precession 

Let us start by defining the vorticity vector, which as usual is given by (in 
relativistic units) 

' e^upu^s (12) 



where the vorticity tensor is given by 

Wa/3 = U[ a .0\ - il[ a Ufl (13) 

and up denotes the four-velocity vector. 

Now, for an observer at rest in the frame of (0), the four- velocity vector 
has components 

Ua = \ '~A'~A~' J (M) 

with 

/V \ 1/2 
A = r-e 2 ? - UWA (15) 

using (0) and 

= e 2/ V 2 sinfl (16) 
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in (0), we easily obtain 

u a = (0,0,0, cj 3 ) (17) 

with 



CO 
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e 



-2/3 o„2/3 4 



{2/3 e e 2/3 - ?Lf< - (uW) 



2r 2 sin $ A 

2f/r 2 2 7 2/3 ([/ r 2 e 2 7 ) jj 2 2 7 

+ -^4— A + A 2 " ~ -45-2^} ( 18 ) 

and for the absolute value of u a we get 

„-2/3-7 /j 

O = (-^ Q ) 1/2 = ^^{2/3 e e^-2e 2 ^-(^r 2 e^) 

2r A v / ^ 

+ 2f/r 2 e 2 ^ + ^ (f/r 2 e 2 ^) - 2/3 u ^[/r 2 e 2 n (19) 
Feeding back into fll9|) and keeping only the leading term, we obtain 

M = -j- ( c ue + 2c„ cot 9) + 0(r" n ) ; n > 1 (20) 

Now, since Q measures the rate of rotation with respect to proper time of 
world lines of points at rest in the frame of (|I|), relative to the local compass 
of inertia, then — Q describes the rotation of the compass of inertia (" the 
gyroscope" ) with respect to reference particles at rest in the frame of (0) (see 
[|nj for detailed discussion on this point). 

Therefore, up to order 1/r, the gyroscope will precess as long as the 
system radiates (c w 7^ 0). Observe that if 

c u6 + 2c n cot 9 = (21) 

then 

F(u) 
sin 9 

which implies 

F{u) = c u = (23) 

in order to insure regularity conditions, mentioned above, in the neighbour- 
hood of the polar axis (sin 9 = 0). Thus the leading term in (|20"D will vanish 
if and only if c u = 0. 
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If the system radiates during an interval of time Am, then the change of 
orientation of the gyroscope, for that period, is given by 

A0 = -VAu (^e^ - U 2 r 2 e 2 ^ V2 (24) 

or, up to terms of order 1/r 

A0 « — (c ue + 2c u cot 6) Au (25) 
2r 

Let us now consider the particular case of a quadrupole radiation in the 
linear approximation. If the quadrupole moment of the source is Q(u), then 
it can be shown (see eqs. (86)-(91) in J7[) that in the linear approximation 

c= l -Q uu sin 2 6 (26) 

and 

2 

- m u = — Q 2 UUU (27) 

Q = ~ sin 26 Q uuu + 0(r~ n ) (28) 

^ 15 sin 26 . , 1/2 ^, _ s ,„„ s 

Q = \/— -m u 1/2 + Or"" 29 

V 8 r 

linking directly the rate of precession to the rate of loss of mass. 

1 



Thus 



or 



4 The gyroscopic preccession of order 

Let us now consider the next order (— ). We easily obtain: 

VI = -—{c u + 2c u cot6») 
2r 

+^ [M e - M{c ue + 2c u cot 6) - cc u9 + 6cc u cot 6 + 2c u c e ] (30) 



7 



Observe that the order — contains, beside the terms involving c u , a term 

not involving news {Me). Let us now assume that initially (before some 
u = Uq =constant) the system is static, in which case 

c u = (31) 

which implies , because of @ 

M e = (32) 

and Q = (actually, in this case Q = at any order) as expected for 
a static field ( for the electrovacuum case however, this may change ||12j| ). 
Then let us suppose that at u = u the system starts to radiate (c u ^ 0) 
until u — Uf, when the news vanish again. For u > uj the system is not 
radiating although (in general) Me ^ implying (see for example @) time 
dependence of metric functions (non-radiative motions 0). 

In the interval u G (uo,Uf) the leading term of the rate of precession of 
the gyroscope is given by (pop. 

For u > Uf there is a precession term of order — describing the effect of 

the tail of the wave on the gyroscope. This in turn provides "observational" 
evidence for the violation of the Huygens's principle, a problem largely dis- 
cussed in the literature (see for example [0],0, [0 and references therein). 

Putting aside the actual technical difficulties in performing such an exper- 
iment, it should be clear that the monitoring of the gyroscope in the interval 
(u < uq, u > Uf) should, in principle, bring out, in a clear-cut way, the 
presence of gravitational radiation. 

Finally, let us consider the particular case of a quadrupole radiation in 
the linear approximation. We obtain for this case (note a misprint in Eq.(87) 
in 0). 

/ 1 K c\ n Of) 1 

n = J^f^f^-^i/a + ±j- 3 Q uuSin2 6) (33) 
V 8 r r z 

Therefore, for u > Uf, the rate of precession is controled by the second 

time derivative of the quadrupole moment (Q uu ). 

5 Scenarios of radiation and estimations 

Let us now present some rough estimates for Q or A0, in different scenar- 
ios. Before doing that, some remarks are in order. 
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Since our intention here is just to provide orders of magnitude of the 
expected effect, we shall restrain ourselves to the quadrupole radiation case. 
Also, although Bondi approach implies axial symmetry, it is not clear that 
such symmetry is present in all examples below. This may be particularly 
true for the case of collisions and bremsstrahlung. In the other cases, specially 
in those of gravitational collapse and supernovae, axial symmetry is not a too 
stringent condition. In the same line of arguments, it should be mentioned 
that in some examples, particularly in binary systems, the signal may last 



for a too long duration. Since the convergence of the series requires (see [10 
for details) 

u<2r (34) 

then, in those examples, the source should be very far from the gyroscope, 
in order to assure the convergence of the series expansion. At any rate we 
should insist on the point that the purpose of this section is not to propose 
specific scenarios for eventual experiments, but just to provide , however 
vague, orders of magnitude of the mentioned effect. 

Now, in the case of quadrupole radiation (in the linear approximation) 
the rate of precession can be related to the rate of loss of mass of the source 
through (^).Then for the change of orientation of the gyroscope A0, during 
an interval of time Au we obtain using 



A0=^^(-m„)^A M (35) 

This last equation will be used to obtain estimations of A0, whenever 
the rate of loss of mass (or the total radiated mass) and the time interval of 
radiation are available. In some examples, when there is not a characteristic 
time scale of the emission we give an estimate of Q, from fl2P|).For simplicity 
the numerical factor and the trigonometric term, in fl35|) and (|29"1), are put 
equal to one. 

In the last three decades, a great deal of work has been done in identify- 



ing possible sources of gravitational radiation (see |14| , and references 
therein). Here we shall present a selection of some of them, keeping in mind 
all reserves mentioned above: 

1. Bynary systems. 
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For the binary system of two neutron stars proposed by Clark [[15 
at 10 Kpc, we obtain A0 pa 2.2 x 10~ 15 . This quantity is obtained 
for an emission time of 1 second, of a binay system emitting mass 
at a rate m u ~ 10 47 Joules/ s. 

For the pulsar 1937+214 (see one obtains fi pa 1.23 x 

10 _23 rad.s _1 . The estimated value used for its angular velocity 
4033. 8rad/s leads to a power of emission around 2 x 10 29 watts, 
taking into account the Landau's formula and considering a dis- 
tance of 2.5Kpc for the source. Identifying the power emitted by 
the system with the loss of mass m u , we obtain the given above 
value of Q . 



• Similar estimations may be done for the pulsar 1913+16 (see |T 
knowing that the power of this pulsar is 6.4 x 10 23 watt. The 
estimates yield, for a distance of 5Kpc, Q pa 1.13 x 10 _26 rad.s _1 . 

2. Gravitational collapse and supernovae 



According to Shapiro fl~5|l , during the first bounce and rebound, 
the efficiency of gravitational radiation never exceeds lO~ 3 M ,then 
for an event at lOKpc with a duration of the order of 1ms (free 
fall time), the maximum obtained A0 is of the order of 3 x 10~ 18 



The model of collapse proposed by Wilson [15|, assumes a radi- 
ating energy (in the form of gravitational radiation) of the order 
of 1O _2 M , during an interval of time of 100M Q (in relativistic 
units). Then for one solar mass , at a distance of lOKpc, one 
obtains A</> pa 6.4 x 10~ 18 . 



In the strongest massive star collapse proposed by Ostriker [[15 
the emission is of the order of 10 _1 solar masses during 10 _1 s, at 
1 Kpc. The resulting A(f> is of the order of 2.9 x 10~ 15 . 

For a model of supernovae proposed by Braginsky and Rudenko 
[|18|1 , (an emission of 10 48 J/s, during 1 ms at 15Mpc) one obtains 
A0 pa 4.2 x 10 -21 . 

The characteristic parameters of a stellar collapse model proposed 
by Rees et al Jl9|, are : an emission of 10 50 J/s during 5 x 10 _4 s. 



In such a collapse there is a continuous frequency espectrum up 
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to — , being t ~ \paGp the characteristic time of the collapse 
t p 

for a final state density p. The variation of energy AE goes like 
1 Q 2 

AE ~ =- where Q is the quadrupole moment along the axis. 

307T T° 

For an event at 5 Kpc involving a star of 6M , the resulting A</> 
is of the order of 7 x 10~ 17 . 

Finally, let us evaluate the emission of gravitational radiation in 
a supernovae event, leading to a neutron star, by estimating the 
total change of the gravitational quadrupole moment of the source. 
Recent estimations of the quadrupole moment of neutron stars 
2(J], point to values of Q of the order of 1.03 x 10 37 kg.m 2 . If 



we assume that the presupernovae massive star has a quadrupole 
moment of the order of the sun ( 3.85 x 10 42 kg.m 2 ) then the 
total change of quadrupole moment is AQ 10 42 Kgm 2 . Taking 
into account fl2~T|) , fl3"5|) , then the obtained A<p, for an event at 
distance of 1 Mpc, during 10 7 years (Kelvin-Helmholtz time scale 
for a star like the sun), is of the order of 2.8 x 10~ 28 . 

Collisions and Bremsstrahlung. 

• The emission, during collision of two neutron stars, given by Wil- 
son [jnjl is of the order of 1O~ 3 M during an interval time of the 
order of 5OM .For an event at lOKpc we obtain A<p « 1.14x 1CT 18 . 

• In the case of two black holes collision, Detweiler [TjJ proposes an 
emission of 1O _3 M , with a collision time of the order of 10M. For 
an event at lOKpc this leads to Acf) « 6.4 x 1CT 25 . 



• For gravitational bremsstrahlung within the galaxy, Ostriker [jT5 
suggest an emission of 1CT 2 solar masses during one second, this 
yields a A0 of the order of 3 x 10~ 16 . 

Other sources. 

• Following a speculative discussion about the evolution of galactic 
nucleus, Ostriker |L5| considers the possibility of destruction of 10 8 



neutron stars at the center, in aproximately 10 6 ' 5 years, emitting 
10~ 2 solar masses in the form of gravitational radiation, an event 
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of this kind, at IMpc of distance would produce a total precession 
of the order of 3 x 1CT 11 . 

6 Conclusions 

We have seen so far that a gyroscope at rest in a Bondi frame will precess 
(up to order 1/r) as long as the system radiates, the rate of precession being 
given by (||). 

Once the radiation stops (vanishing news) the gyroscope will continue to 
precess with a rate of rotation given by the second term of (|30|) with c u = . 

As can be seen from estimations done above, excluding the last example, 
and the example based on the evaluation of the change in the quadrupole 
moment, the most realistic scenarios point to A0 's ranging between 10~ 15 
and 10~ 19 . We ignore how far are we , with the present technology, to the 
accuracy required for that kind of measurement. Nevertheless, we want to 
stress that it has been our main purpose here, just to bring out such effects 
in the context of Bondi formalism. 

We would like to conclude with the following comment: observe that all 
along our discussion we have not made reference to specific bandwiths. This is 
so because all quantities in the Bondi approach, are not defined with respect 
to any specific frecuency. In this sense eq. (|TT1) (and all related quantities), 
have to be considered as integrated over all frecuencies. 
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